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ABSTRACT 
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1 INTRODUCTION 


From Newton’s law of universal gravitation to Einstein’s general relativity, our under¬ 
standing of gravity has developed profoundly. However, efforts to incorporate quan¬ 
tum mechanics into gravity have been unsuccessful, at least partly due to the fact that 
quantum gravitational effects are unobservably weak at low energies jlj. This unsat¬ 
isfactory situation has been transformed by the accumulation of large cosmological 
data sets on primordial perturbations which are predicted by the theory of inflation 
to arise from the gravitational response to quantum fluctuations of matter and from 
quantum fluctuations in gravitational radiation [2], [3], 01 [5], El 0 IH1 El EES EEE1 H2] . 
These quantum gravitational effects are hugely enhanced during primordial inflation 
and then fossilize so that they can survive to much later times [13]. 

Electromagnetism provided humanity’s first and simplest example of a relativistic 
unified field theory. It is therefore natural to wonder how electromagnetism affects 
quantum gravity during inflation. Of course, Einstein plus Maxwell is not perturba- 
tively renormalizable, even at one-loop order mm, and that does mean we must 
accept some limitations on what can be studied with it. We use Einstein plus Maxwell 
in the standard sense of effective field theory, absorbing ultraviolet divergences order- 
by-order with Bogoliubov, Parsiuk, Hepp ,and Zimmermann (BPHZ) counterterms 
[16]. Once this is done, the renormalized results from loops of massless gravitons 
and photons engender nonlocal and ultraviolet finite contributions to the effective 
action which are unique predictions of the theory that cannot be changed by its still 
unknown ultraviolet completion |T7, IS]. A famous example is the computation by 
Bjerrum-Bohr of the one graviton loop correction to the Coulomb potential on a flat 
space background |I9]. We work in the same theoretical context but on a de Sit¬ 
ter background, and with the slightly different goal of inferring the one photon loop 
corrections to the Newtonian potential. 

One way of studying quantum weak-field corrections to a classical theory is to 
first compute the corresponding one-particlc-irreducible (1PI) two-point function. 
This is then used to quantum correct the linearized effective field equations. For 
our problem we must compute the one photon loop contribution to the graviton self¬ 
energy —i [ /U 'E P<T ] (x; x'). We can then infer quantum corrections to graviton mode 
functions, or to the gravitational response to sources, from the linearized effective 
field equation for the graviton field h pv (x). 


V^Kh pa (x) 


(x-x')Kh pa (x') = TZ(x) = 


5S, 


matter 


5h pu (x) 


( 1 ) 


h =0 


Here, 1~)^p ct is the Lichnerowicz operator in a de Sitter background, k 2 = 167tG is 
the loop counting parameter of quantum gravity, and T^{x) is the linearized stress 
tensor density of whatever source is desired. In our case this source is that of a 
static point mass. The same framework (with zero source) has been employed to 
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investigate corrections to graviton mode functions from massless, minimally coupled 
(MMC) scalars [2D, [23]. There have also been studies of the effect of a loop of 
massless, conformally invariant scalars on the graviton mode function [22, (22SJ and on 
the background geometry [2Tj . 

This paper contains six sections, of which the first is this introduction. Section [2] 
gives those of the Feynman rules for Maxwell plus Einstein which are needed for our 
computation, as well as some facts about the background geometry and the way we 
represent the tensor structure of the graviton self-energy. In section |3l we derive the 
primitive, dimensionally regulated result for the one photon loop contribution to the 
graviton self-energy. Section [4] introduces appropriate BPHZ counterterms for this 
theory and gives the result of renormalization. In section [5[ we convert our in-out 
graviton self-energy to the retarded one of the Schwinger-Keldysh formalism, and we 
use it to solve ([I]) for a stationary point mass. In section [6] we propose a physical 
interpretation of our result. 

2 FEYNMAN RULES 

2.1 Preliminary clarifications 

This work is based on quantum gravitational perturbation theory in de Sitter space 
which is known to be the unique maximally symmetric solution to the Einstein held 
equation with a positive cosmological constant A. But we will not work on the full 
D-dimensional de Sitter manifold rather on a submanifold known as a “cosmological 
patch” which is homogeneous, isotropic, and spatially hat. The invariant line element 
for this geometry is 

ds 2 = —dt 2 + a 2 {t)dx ■ dx , (2) 

where aft ) is the scale factor measuring the expansion of the Universe. The Hubble 
parameter is defined as H{t ) = d{t)/a(t). It is convenient to use open conformal 
coordinates U' = (?/, x l ) with 

— oo < rj < 0 , — oo < x 1 < +oo , for i — 1, 2, • • ■ , D — 1 , (3) 

where the conformal time coordinate 7) is related to the normal time coordinate 
t through the relation, dt = adi). The metric in open conformal coordinates is 
conformally hat, 


ds 2 = a 2 (;if) {—dr) 2 + dx ■ dx) = a 2 fr))ri^dx^dx u . (4) 

The Hubble parameter H is constant in de Sitter spacetime, and the scale factor in 
terms of the conformal time coordinate is a(i]) = — -A-. We then define the graviton 
held h^fx) by conformally transforming the full metric g lw {x) and then subtracting 
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Figure 1: Photon contributions to the one-loop graviton self-energy. Photon propa¬ 
gators are wavy and graviton propagators are curly. 

off the background, 


9^(x) 


a 2 (v) Vttv + Kh 


a 2 (r])g^(x) , 


(5) 


where rj /lL/ is the D-dimensional Minkowski metric with a spacelikc signature and h^ u 
is the graviton field whose indices are raised and lowered with the Minkowski metric. 
And we will make great use of the de Sitter invariant length function y(x\ x r ) while 
representing our propagators on de Sitter spacetime, 


y(x;x') 


a(g)a(g')H 2 ||x — x !|| 2 — (|?/ — rf[ — i8)“ 


( 6 ) 


2.2 Primitive diagrams 

All possible propagators and interaction vertices are derived from the primitive La- 
grangian density £, 


C 


£gR + £eM + ^GRfix + ^EMfix + ^BPHZ • 


(7) 


Here, the gravitational Lagrangian density of low energy effective field theory is 


1 r 


£gr — ~ 


K 


R — (D — 2)A y/^g , 


( 8 ) 


where the factor of (D — 2) multiplying A makes the pure gravitational held equation 
take the form of R^ u = Rg fW in any dimension. And we define the cosmological 
constant as A = (D — 1 )H 2 . The Lagrangian density for electromagnetism takes the 
familiar form of Maxwell’s theory, 

£em = —F^F^g^s^Tg , (9) 


• Propagators 

To get the graviton propagator, we expand the Lagrangian density ([8]) with the full 
metric (J5]) and extract a presumably irrelevant surface term. Then, the quadratic 
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part takes the form 


£gr — Surface Terms = 


D 


- l) Ha D - l s/^gP^ u K 0 d,h pa 


To \/' ~9 9 9^ 9 ^ ^ 2 ihapdfjhi/a- ^ d p h fy j c) iy h (T p ^ d p h (y '^() (J h pu ^dph ap d v hp 0 


We fix the gauge by adding 
1 


£cRfix — 


D - 2 < V F»F V , F p = rf a (dji fip - l -d p h pa + (D - 2)Ha5° a h^ . (11) 


This gauge fixing term (TlTj) is nice because it allow us to express the graviton prop¬ 
agator as a sum of scalar propagators multiplied by constant tensor factors, 


i [ab/AoJ (x; x') = xfA 7 (x;x'). 

I=A,B,C 


( 12 ) 


The various tensor factors are 
= 2 %</}*)» 


rpA 
V-v-*- pa 


(.0-3) 


VpvVpa i 


rpB 

pa 

rpC 
v-v 1 pa 


~^(pVu)(p S a) 1 
2 

(£>-2)(£>-3) L 


(£>-3)^ + ^ (D-3)^X + V 


(13) 

(14) 

(15) 


where r) lxl ' = r) tlv + SqSq is the spatial part of the Minkowski metric. Each of those 
scalar propagators obey the associated propagator equation, 

UjiAi(x] x') — iS D (x — x') , I = A,B,C , (16) 

where we define three scalar kinetic operators as follows: 

V A = dp (y/^srdv) = □ , (17) 

and 

V b = V a -{D- 2 )H 2 ^ , V C = V A - 2{D - 3)H 2 ^9 ■ (18) 


.( 10 ) 
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Note that iA A (xyx') is just the propagator of a MMC scalar. Thanks to previous 
work [2S], all of the scalar propagators have been worked out, 


iA A (x]x') = iA c f(x]x') 
H d ~ 2 T(D - 1) 


(4vr)# r(f) 


D 


d r 2 (f) /4^f-2 

AT(D-l)[ yj 


( 7T 


7rcot |^—DJ + 111 ( 00 ') 


h d ~ 2 

1 T(n+D-1) /y\ n 1 T(n+f + l)/ 

1 (4tt)t 

n T(n+f) \4/ n-f + 2 I>+2) ' 

*Ab(x; x') = B(y ) = iA cf (x;x') 

jjD-2 00 

fr(n + L>-2) r(n + f) /?/\"-#+ 2 l 

( 4 +>^ 5 

[ T(n+f) V4/ T(n + 2) V4/ J 


+-f + 2 


iA c (x-,x') = C(y ) = iA cf (x;x') 


H 


D—2 


(47T 


X + ( n + !) 


T(?r + D — 3) fy 


n =0 


r (n + f) V4 


n 


D 

2+ 3 


r(n + f - 1) ry 
r(n + 2) U 


(19) 


( 20 ) 


( 21 ) 


where iA cf (x; x') is the de Sitter invariant propagator of a conformally coupled scalar, 

H d ~ 2 fD \ /4\^ _1 

tArf(liI ^_ r (__ i )(_) . (22) 

Note that the inhnite series in all of these scalar propagators vanish in D = 4 dimen¬ 
sions, which means they are only needed when multiplied by some divergent pieces. 

Unlike gravitons, photons show no physical breaking of de Sitter invariance. So, 
we conlcl in principle employ an exact de Sitter invariant gauge to make the photon 
propagator manifestly de Sitter invariant [27]. However, a noncovariant gauge fixing 
£embx 123 !2Bj results in a much simpler form of the photon propagator, 

i [ m Aj,] (x; x') = rj^ x aa'iA B (x] x') - x aa'iA c (x ; x') . (23) 

Here and henceforth we adopt a shorthand notation for the conformal time depen¬ 
dence of the scale factor, 

a = a{rf) , a' = < 2 ( 7 /) . (24) 


• Vertices 

Note that the action of electromagnetism, S'em = jd D xC B Mi contains interactions 
between two photons and any number of gravitons. So we obtain the three-point 
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and four-point interaction vertices by functionally differentiating S'em once and twice 
with respect to the graviton field, and then setting the graviton field to zero, 


^EM 

Sh^x) 


5 2 S( 


EM 


(^)^^pcr (•£ ) 


h=0 


h=0 


Ka D - 4 V aK ^d K A a (x)d x A p (x) , (25) 

K 2 a D ~ 4 U a ^ xpvprj d K A a {x)d x A p (x)5 D {x - x') . (26) 


The tensor factors for the three-point and four-point vertices are [29] 

yOK^XpU _ ^ ^X]k 4 r jV)[°‘ r jK][P r jM( l/ ^ 
jjanpXpvprr _ ,^a[^A]re 

-)-r7 /i!/ 77 p )[ a 77 K l[^7/ A ]( <T + rj pa r)^ a r] K ^r] X ^ 1 ' 

_| _^a{p^v){p^a)[^X\K, _|_ n(p v)(p a)[X p]a 

_)_^«(m^^)[/3^A](p^(t)k n(p v)[X 0](p <r)a 

_|_ ^«[A £]0 ? ^)(p^v)a _ 


(27) 


(28) 


• Diagram Expressions 

Now, using the Feynman rules stated in section [21 we can write down the formal 
expressions for the first two diagrams in Fig. [I] The leftmost diagram, also known 
as the “three-point function”, is constructed from two three-point vertices, 

-• ['“'£&] (z; X 1 ) = i (-in) a 0 - 1 [„A„] (x; x') 

x (-in) a’ D - 4 V^ pa ded^ [ 7 A 5 ] (x; x') . (29) 

The middle diagram, also known as the “four-point function”, is constructed from a 
single four-point vertex, 

- i [^E(£ t ] (x; x') = (-in 2 ) a D ~ A U aK ^ upa d K d' x i [«A^] (x; x')5 D (x - x') . (30) 

The rightmost diagram denotes a series of counterterms which will be explained in 
section 14.11 


2.3 Structure functions 

The graviton self-energy possesses 100 components in D = 4 dimensions. We could 
report on all of them, but the symmetries will relate these components such that 
they can be expressed using a few structure functions. This idea is actually inspired 


6 








by experience in flat space. For example, the Poincare invariance makes the vacuum 
polarization of photons take the form of a single structure function II((x — x ') 2 ), 


And the graviton self-energy takes the form 



U pu U pa F 0 ((x-x' 


^YIv{pyy?)v 


n^n^\ 


F 2 ({x-x') 2 ) , (32) 


where we dehne the projection operator as II= d p d v — ij^d 2 . We usually call F 0 
and F 2 as the spin zero and spin two part. The de Sitter background for inflation 
is spatially homogeneous and isotropic. Then, if these are the only isometries, we 
take it to be manifest in our representation for the graviton self-energy that it turns 
out that we shall need four structure functions. And guided by the previous studies 
[301 |3I] about a noncovariant representation of the vacuum polarization on de Sitter 
spacetime, 

i [^hP] (x; x') = (rTrTd p & a F(x] x') + {fj^fj pa -fj^fj vp ) d p & a G{x\ x') , (33) 

which makes interpretation of the physics much easier. Here, we will adopt a similar 
noncovariant representation of the graviton self-energy [32j , 


-i [^S pfT ] (x; x') = F pu (x) x F^^x') F 0 (x]x') 


FG^ix) x Q pa {x') G^x-x') 


_j_ jrpvpa 


F 2 (x-,x') 


_j_ gpvpa 


G 2 (x-,x') . (34) 


where two scalar projectors, F pu and G ,w , are 

F pv = d p d v + 2 (D - 1 )aHS i 0 lx d u) + (D - 2)(D - 1 )a 2 H 2 8 p 8" 

-rf v [d 2 + (D- l)aHd 0 + (D - 1 )a 2 H 2 ] , (35) 

G ,w = + 2 (D - 2 )aH5^d v) + (D - 2)(D - l)a 2 H 2 8$6% 

-1f u [V 2 + (D - 2)aHd 0 + (D - 2 )a 2 H 2 ] , (36) 

and where = d p + And two tensor projectors, F pvpa and Q pvpa , are 

T pvpa = Q/ i/ WxV(^)x^V\ 7 V, (37) 

G pvpa = VWx^fi'JxfyV^, (38) 

where we dehne the second order differential operator C a g y s pu by expanding the Weyl 
tensor of the conformally transformed metric to the linear order of the graviton held, 

9pv = 9pv + G a gry$ = x hvh^ LU A G(k h ) . (39) 


7 












Explicit expressions for F^ ppa and Q pup(T can be found in the Appendix of [32] . 

ft is worth mentioning that the de Sitter representation (1311) must be consistent 
with its corresponding flat space limit (I52|) . This means that all terms in Gq{x\x') 
and G 2 {x-,x') must contain at least one factor of H 2 , so as to vanish in the limit 
that H vanishes with the comoving time t = In {a)/FI held fixed. Hence, the lead¬ 
ing divergences are contained in Fq(x\x') and F 2 (x] x r ). Gq{x\x') and G 2 (x]x') are 
correspondingly less divergent. 

Now this leads to the question, how do we find the structure functions? In general, 
suppose that after substituting the photon propagators in our formal expressions given 
by Eqs. (129|) and (130|) . doing index contractions and acting derivatives, we obtain a 
primitive result for — i ^E pcr ] (x; x'). Since this result can be expressed in the form 
suggested by Eq. (1341) . as we explained previously, we can extract out those four 
structure functions by looking at particular components. The procedure will be as 
follows. 

• First, we take the trace on one index group of — i [^E^], say p and a. This 

will cause the spin two part to drop out. Then we shall look at the components 
—i [°*E p,T ] x r] pa and —i x r] pa (j ^ k ) to obtain two linearly independent 

equations to solve for F 0 (x',x') and Gq(x;x > ). 

• Next, we look at the components —and —i [t fc E 0 *], and always assume 
that i 7 ^ j 7 ^ k i. This will also allow us to derive two linearly independent 
equations to solve for F 2 (x]x') and G 2 {x^x'). 

More details about how this procedure works will be given when we carry out the 
explicit calculation in section 13.21 and 13.31 

3 ONE-LOOP GRAVITON SELF-ENERGY 

3.1 Primitive results of three-point and four-point functions 

In this section, we will substitute our vertex tensor factors (127|) and ([251) into the 
formal expressions (129|) and (130|) and then perform naive contractions. Before diving 
into the tensor algebra, it is worth noticing the fact that in both of the diagram 
expressions (1291130]) the doubly differentiated photon propagators inherit an antisym¬ 
metry from the vertex tensor factors through index contractions. So, we will have a 
cleaner, simpler analysis if expressing the computation in terms of a four-index object 
which is the correlator of two field strengths, 

i[ a pA^ s ](x-,x') = (fi|F aj8 (x)F 7 i(a/)|fi) = ^d [a d' [[ry i A 5]] ] (x; x') . (40) 

Here, the single and double square brackets indicate antisymmetrization on the index 
pairs a -H- (3 and 7 Pi5, respectively. 
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Now let us evaluate the field strength correlator (1501) . Starting with our photon 
propagator expression (1251) . we want to first express the tensor factors in terms of the 
derivatives of y(x] x') and u{x\x') = 111 ( 00 '), 


««V = 


2 H 2 
1 


= 2Hn 


-d,,dly + d p y x d' v u + d p u x d' v y + (2 - y)d p u x d' v u , (41) 

(42) 


2 dfjU x & u u 


And there is a very useful relation between B(y ) and C(y), 


Cto) = 5(2-#)fl(») + ^ 

Substituting expressions (15T1) - (143|) into 


where k ^^T ( D-l) 


(43) 


(4w)f r(f) ' 
allows one to express our propagator as 


(x;a/) = 


2H 2 


dpdljy- d p yd' v u - d p ud' v y 


k 


(D 


3 } H 2\ d » ud 'v U 


(44) 


The next step is to act two derivatives on (1551) and antisymmetrize them to form 
the field strength correlator (j50l) . Notice the fact that d p u = Ha5 p is independent 
of and that d' u u = Ha'5® is independent of x^. It follows that only the de Sitter 
invariant first term of expression (155)) contributes to the field strength correlator, 


i [ a/J A 75 ] (:r; x') = (. D [a D[ b B(y )) (D^D^yj , (45) 


where D a and D' are the covariant derivatives associated with the de Sitter back¬ 
ground metrics p (:r) = a 2 7] pi/ and g pc (x') = a' 2 y pcT1 respectively. It also turns out to 
be useful to work out the coincidence limit of the field strength correlator. Starting 
with (155|) and using the fact that 


lim (D a D'B(y)) =-2H 2 a ,2 B'(0)r] ary where £'(0) = limf^-SQ/)) . (46) 

x—>x' \ ' J y —>0 V ay / 

Hence, the coincident field strength correlator is 

lim i [ Q/3 A 7<5 ] (x\ x') = —8 H 2 a ,A B'(0)t] a ^r]s}i3 ■ (47) 

x—>x' 

The final identity we would like to derive for future convenience is based on the fact 
there are, of course, local delta function contributions hidden in the expression (1551) 
because of the doubly differentiated scalar propagator. And one can show that 

D p D' u B(y) = {D p D' u B(y)\ + %8° v i6 D (x - x?) , (48) 

^ J naive CL 
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where {.. .} naive denotes whatever we got by naively acting derivatives. Then, the 
result (USD allows us to segregate out the local delta function contributions in the 
field strength correlator, 



Now we are all set to compute the primitive results of the first two diagrams in Fig. |T| 
by substitution of (127|) and (1281) in (1291) and (130|) and naive index contractions. First, 
it is easy to see that in terms of the field strength correlator (ITUD the four-point and 
three-point functions can be rewritten as 



Ji/’i [“"A 7< ] (z; x') x i P'Atj] (x;x') +i [<“'‘A 7(p ] (x-,x')xi [VA/ 1 ] (z; x') 


(51) 


Note that while getting (j50l) and (T5TT) we raise or lower indices only with the Minkowski 
metric, which just plays with the relative minus sign between the temporal and spatial 
components. 

• Local contributions 

Note that all terms in (1501) are proportional to the delta function. So, applying 
the coincident identity of the field strength correlator (1471) . we immediately obtain 
the contributions from our four-point functions which are purely local, 



(52) 
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Obviously, the expression (j5Tj) will contain both local and nonlocal contributions, and 
we first deal with the local terms. So, applying the identity (jUJj) and ignoring all of 
the naive terms for the time being, we can also easily obtain the local contributions 
from our three-point function, 


= ~^H 2 a' D B'(0)l(D - 5)rTrT + W’" 

+2 (D - 4) ij'“'W+ 2 ( b - 4) (D - 4) \ti D (x - x') . (53) 


Adding (T52|) and (J53J) together simply gives the entire local contributions from our 
three-point and fonr-point functions, 


-rsn„e>; *') = (D ~ 4)« 2 ifY D B'(0) -ID- OW" 


-Ud - 5)r- 2rTK% ~ 2i ftiZK + 4$^$ \iS D (x - x') . (54) 


From our explicit expression of B(y) (j20|l . we know that B'{ 0) is finite in D = 4 
dimensions, 

B'(0) = . (55) 

4(4ir)fr(f + 1) ' 

So, we can simply set D = 4 and conclude that the net contribution of (l54]i vanishes. 


• Nonlocal contributions 


Now let us focus on the nonlocal contributions which are nothing but all of the naive 
terms, that produce no delta functions, in our 3pt function. Please note that (|5Tjl is 
not written in a manifestly covariant form because we have raised indices using the 
Minkowski inverse metric instead of the inverse background metric. That is merely 
for the purpose of naive index contraction, but we can easily recover the covariance by 
putting back the appropriate factors of a or a'. Also, note that the tensor structure 
of the naive terms in our held strength correlator (l45]i can be expressed in terms of 
the covariant basis furnished by derivatives of y(x;x'), 


i [ ap A 7<5 ] (x; x') 


2 r 

JP 


(D [a y) {DflD' b y) ( D' 5] y ) B"{y) 

-{D a D' h y) {D' s] D p y) B\y) 


(56) 
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And one can show there are some very useful contraction identities of these basis 
tensors [28] , 


(D„y) (D„y) = H 2 (4y - y 2 ) = g^(x') (C» (D» , (57) 


r l 'hri (D„y) (D v D' f y) = H 2 (2 - y ) (l/,,</) , (58) 

rV) (D'„y) (D r D' a y) = H 2 ( 2 - y ) (D,y) , (59) 

9 , “'W (r>,C» (A^i,) = 4i/X(i') - ff 2 (D' f y) (D'„y) , (60) 

rMMlDA) = iH l g li .(x)-H 2 (D,y)(D v y) . (61) 


We remind the reader that = a^ 2 g pi/ and g pa = a'~ 2 ij p(7 are the inverse background 
metrics on the de Sitter background. Next step is to substitute our result ()56ll for the 
naive terms in the expression (ED), and then to make use of the contraction identities 
(l57)b([6Tf . ft is straightforward to express all of the naive terms in our three-point 
function as a linear combination of the de Sitter invariant basis tensors, 


-i[ pu Z prT ](x;x') = (aa') D+2 ^(D p D'( p y) (D'^D u y)a(y) 


+H 2 


r u (x) {D'^y) (D^y) + g pa (x') (D^y) (D^y)] e{y ) l . (62) 


And those five scalar functions in (16211 are given by, 


a 


(</) = -f (4;/ - y 2 Y [s"(</)] 


2 (D - 2) - - (4y - y 2 ) 


[B'(</)f 


--(2 - y) (4y - y 2 ) [B'(y)B"(y)} \ , (63) 


m = k2 1 4 < 2 - V) (4</ - y 2 ) [B"(y)} 2 - (2 - y) [B'(j/)] : 


2 (D - 4) + (iy - y 2 ) 


[B\y)B'\y)\ 1 , (64) 


y(y) = k 2 


\(D -2)-\(4y- y 2 ) 


[B"(y)} 2 ~^[B\y)\- 


--(2 -y)[B'(y)B" (y)] }- , (65) 
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<5(3/) = k. 2 {- „(D - 5)(4y - yJ[B'\y)] 2 - 


( D /_9£> + l 6 )__(£(_5)(4 y-y 2 ) 


{B'(y)}‘ 


- 2<-0 - 5)(2 - y) (4y - y 2 ) [S'(y)B"(y)] f , ( 66 ) 


iv) = k 2 \ j (D - 4) (4y - y 2 ) [B"(y )] 2 - (D - 4) [fl'(y)] : 


+(B — 4)(2 — y)[B'(y)B"(y)]| .(67) 

Since our local contributions (154|) vanishes at one-loop order, the nonlocal terms (l62j) 
complete the one-loop photon contributions to graviton self-energy. 

3.2 Spin zero structure functions 

In this section, we will follow the procedure introduced in section 12.31 to extract 
two spin zero structure functions F 0 (x]x') and Gq(x]x'). We start with the trace of 
expression (E2D, 

-i pa } (x; x') x r lpa = a 2 H 2 (aa') D ^H 2 r^x) [ 4 a(y) + D8(y) + ( 4 y - y 2 ) e(y) 

+ (D^y) (D v y) -a(y) + (2-y)(3(y) + {Ay-y 2 )^{y) + De(y) J. ( 68 ) 

Now, recall that our strategy is to look at the components —*[°*£ pf7 ] x rj pa and 
—i x y pa (j 7 ^ k ) to infer equations that F 0 and Go satisfy. One can see that 

the time derivative on y(x]x') gives 


d 0 y = aH \ y- 2 + 2—J , (69) 

which leads to the components we need of (D^y) ( D u y ), 

(D°y) ( D l y ) = -a -4 (d 0 y) (, d®) = -a~ 4 x aH (^y - 2 + 2^j (, d t y ) , (70) 
(D J y) ( D k y ) = a -4 x (dfy) (d k y) . (71) 

As explained in work [32], now we need to recognize two important scalar sources 
Si(x]x') and S 2 (x-,x'), by writing 

— i [°*E pCT ] [x] x') x r] pa = a D ~ 2 a' D x diSi(x\x') , 

—i p fc S pcr ] (x; x') x rj pa = a D ~ 2 a' D x djd k S 2 (x\x') . 
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(72) 

(73) 








If we define the indefinite integral of a function f(y) as I[f(y)\ = f y dzf(z), and 
another scalar function F(y) as 


F(y) = I 2 


-a(y) + (2 - y)0(y) + (4y - 9 2 ) 7 ( 9 ) + £>£( 9 ) , 


(74) 


one can see that 

-i[ 0i ^](x-x') xr, pa 

-i[ jk ^]{x-,x')xri pa 
It follows that 


a D 2 a ,D dA—aH 3 I 


0 V ~ 2 )F"(y) 


-2 a'H 3 F' 


a D - 2 a' D d 3 d k {H 2 F(y)} . 


(l/)} X75) 

(76) 


Sifa a/) = H 2 (-d 0 + aH)F(y) , (77) 

S 2 (x-x') = H 2 F(y) . (78) 

Then, substituting (177)1 and (1781) into (56) and (57) of [32] allows us to express the 
spin zero structure functions as 


F 0 {x]x') 

G 0 (x]x') 


-MIT! (-*—) F(y) 
(D - l) \U + DH 2 ) ’ 

0 , 


(79) 

(80) 


where □' denotes the invariant scalar d’Alembertian evaluated at i.e., □' = 
g pa (x')D' p D' c = a'~ D d p (a' D ~ 2 r] pa dP) ■ Now, let us first work out what explicitly the 
scalar source F(y) is. Substituting (I63|) - (l65|) and (167)) into (174|) gives 


F(y) = P(D-2)(D-4)I 2 


(■ B'f + (2 - y) (B'B") + i (49 - y 2 ) (B"f 


(81) 


Then, substituting the explicit expression of B(y) (120]) leads to the following expan¬ 
sion of F(y): 


F(y) 


(. D - 4) 2 K 2 H 2D ~ l T 2 (f) 
II,(D - 1) (4ir D ) 


(D-2) 



(D- 4) 



. (82) 


where we have neglected terms that are integrable and vanish in D = 4 dimensions. 
Similar terms will also be regarded as irrelevant in subsequent expansions. 

The next step is, of course, to recover the structure function F 0 (x;x') from ex¬ 
pression (1791) by employing the Green’s function for (□' + DH 2 ). But note that 
electromagnetism is conformally invariant in D = 4 dimensions, which means that 
our graviton self-energy shall be traceless in those dimensions; therefore, the spin zero 
structure function F 0 (x ; x r ) is expected to contain a factor of (D — 4), as we see in our 
result of source term F(y) (1821) . However, there happens to be a factor of (D — 4 ) 2 
in F(y) which causes absolute zero one-loop contributions to F 0 (x\ x') since the most 
divergence we can have at this order is (D — 4)~ l . To sum up, both of the spin zero 
structure functions, F 0 (x;x') and G 0 (x]x'), vanish at one-loop order. 
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3.3 Spin two structure functions 


In this section, we intend to extract the spin two structure functions F 2 (x\ x') and 
G 2 (x]x') from the components — i (x; x') and — i p fc E 0 *] (x\ x'), with i ^ j ^ 

/■• / / m- 

(83) 

(84) 


— i[ 0l Y, jk ](x]x') = ( aa') D 2 djdjd k S 3 (x] x') , 
—i (a;; x') = {aa') D ~ 2 didjd k S^{x] x') . 

This means that we will need the following identities: 


{D°y) (D'y) = - a ' 4 x aH [y - 2 + 2-J (d iV ) , 

(D'°y) ( D H y ) = a '" 4 x a' H [y - 2 + 2 ^) (<^) , 

(D°D ,i y) = ( aa ')” 2 x a// (3#) , (tW?/) = -( aa ')~ 2 x a'# (<%?/) 

Applying these identities to the desired components of gives 


(85) 

( 86 ) 
(87) 


-i[“ T. jt ](x-,x') = (aa') D ^a l d J d t \--HaI 3 [lS(y)\-(d 0 -3Ha)I i [-i(y) ^,(88) 


-i P*E“] (x; x') = (aa') D - 2 d, d,d k \-Ha'I 3 [/?(;,)] + (SJ, -3/fa') /* [ 7 (») 


Then, comparing the results 
other two scalar sources, 


and ()89l) with expressions 


and 


1 


(89) 
gives an- 

(90) 

(91) 


Substituting expressions (l90jl and (l9Tf into (66) and (67) of [32] and taking into 
account that in our particular case both F 0 (x] x r ) and G 0 (x',x') vanish, we find that 
F 2 (x]x') and G 2 (x]x r ) are determined by the following two independent equations: 


S 3 (x;x') = --Hal 3 /3(y) - (<9 0 » 3 Ha) I 4 ~/(y) 
S,(x-x’) = ^Ha'I 3 [l3(y)]+(d' 0 -3Ha')I^( y ) 


(do - do) 


(D- 3) _ 

(D-l) 2 (D-l)(D-2) 


(D- 3) 


:G, 


(aa') D ~ 2 I 4 7 (y) 


= (a - a')H(aa') D ~ 2 x { --I 3 [(3(y)\ + (D + 1 )I 4 [i(y) 


(92) 


(d'o + do 


(D-3)D 


-F 2 + 


(0-3) 


G 2T(on ) 


l\D- 2 J -4 


(D — 1)(D — 2) z (D-l)(D-2) 2 
= (a + a')H{aa') D ~ 2 x (-i / 3 f/3(j/)l + (D + l )/ 4 [ 7 ( 1 /) 


7 ( 2 /) 


(93) 
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To solve equation (j92j) . consider acting {d' Q — d 0 ) on ( aa') D 2 f{y), 


(<% - d 0 ){aa') D ~ 2 f{y) = (a - a')H{aa') D ~ 2 (4 - y)f{y) — (D — 2 )f{y) 
If we let 


(aa') D -y(y) = 


(0-3) 


;G 2 — {aa 1 ) 


l\D -2 tA 


{D- 1) z (D — 1){D — 2) 2 
then (l92]l becomes a first order linear differential equation, 

1 


7 (y) 


f'(v) - -F—lf(y) = 1 


~i 3 m +(o+ 1 )/* 7 (») 


(4 — 1 /) (4 - y) 

which allows us to invert (dg — <%) in (1921) and obtain an indefinite integral, 


(94) 


(95) 


(96) 


(P-3) 

(P-1) 


F 2 (x; x') 


(P-3) 


(aa‘ 


t\D -2 


( 4 -y) 


D—2 


(D — 1){D — 2) 2 
1 


G 2 {x; x') — {aa') 


l\D-2 tA 


7 (y) 


(4 -y) 


D -3 


m + <p + i)p 7(») 


+ M , (97) 


where the integration constant K\ is chosen so that there is no singularity at the 
antipodal point y — 4. Similarly, to solve Eq. (j93j) . consider acting {d' Q + do) on 
{aa') D ~ 2 f (y), 

(d'g + d 0 ){aa') D - 2 f{y ) = (a + a')H{aa') D - 2 [yf'{y) + {D — 2 )/(y)] , (98) 

which allows us to invert {d' 0 + <9q) in (1931) and obtain another indehnite integral, 


(D 3 )D . .. {D 3) , / „/\D-2 t a 


(D-1)(D- 2) 


F 2 (x;x') + 


(aa 


r\D -2 


y 


D—2 


{D -1){D -2) 2 
1 


G 2 (x; a/) + (aa') J 


,D—3 


-I 3 0(y) +{D + 1 )/ 4 7 (y) 


7(2/) 

+ /i 2 , (99) 


where the integration constant K 2 is chosen so as to keep G 2 {x\x') less divergent 
than F 2 (x]x'). Equations (1971) and (j99|) simply imply that 


F 2 (x;x') = 


G 2 (x-,x') = 


{D — 2) J {aa') D 2 „ (aa' 


l\D-2 


2{D-3)\{A-y) D ~ 2 
{D — 1){D — 2) 2 f (D-3) 


—Pfe) + —5—<?(!l) i . 


( 100 ) 


(P-3) 


(P-1) 


F 2 {x-,x') 


(aa 


J\D -2 r 4 


7(2/) 


aa') 13 2 ... 

J -P(y) . (101) 


(4-2/) 


D—2 
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where we have defined two functions of y(x] x'), 

1 


P(y) = I 


Q(y) = 1 


( 4 -y) 


D -3 


--J 3 0{y) + (D + l)/ 4 7 (j,) 


y 


D -3 


--I 3 

2 


P(y) + (a> + !) /4 7(2/) 


Substituting the explicit expression of B(y) (l20j) into (164)) and 
various indefinite integrals, we obtain the following expansions: 


i 3 [m 

r'bt v) 


2 u2D-4p2 / ' D 


k 2 H 


(3 D ~ 8) 


(47t) D 

k 2 H 2D-4 T 2 ( D ^ 


(4vr; 


D 


200-2)00-1) Vy 

0> / 4 

'800-1)00 + 1) U 


+ Ad , (102) 

+ K ‘2 . (103) 

and working out 

(104) 


D—2 


D—2 


(105) 


As always, the neglected terms in these expansions are integrable and vanish in D = 4 
dimensions. Thanks to the factor of (D — 4) 2 in the following result, which is a little 
surprising, 

1 


--/ 3 [/%)j +(D + l)I^(y) 

2 u2D—4p2 ( D 


k 2 H 


(■D - 4) 2 


(47 t) d 1 8(D-2)(D- 1) \y 

our analysis becomes a lot easier. As a result, we find 

P(y) = K 1 , Q(y) = K 2 . 


D—2 


(106) 


(107) 


Note that we shall set K\ = 0 to remove the singularity at the antipodal point y — 4. 
Then, it follows that 


F 2 (x-,x') = 


G 2 {x-,x') = 


(D — 2) I (aa!) 


l\D—2 


2(A> — 3) I y 


D—2 


-K, 


2 ( , 


(108) 


(D — 1)(D - 2) 2 (D — 2) (aa') 


l\D-2 


(D-3) 


2(B-1) y 


D—2 


-IU. 


K 2 (H 2 aa') D -' 2 T 2 (A) 


D 


(4vr; 


D 


8(A> — 1)(A> + 1) \y 


D—2' 


(109) 


Note that (jl08j) and (11091) suggest that both F 2 and G 2 are proportional to l/y D ~ 2 . 
Recall that G 2 should be less singular than F 2 : therefore, the integration constant 
K 2 can only be chosen so that G 2 vanishes at this order, 

2„ D—2 _ p2 ( D\ ( / 4 \ D—2 ' 


(aa’) D 2 _ K 2 (H 2 aci‘ 

y D ~ 2 2 (47t 


D 


\D 


4(A) — 2)(A) + 1) \y 


( 110 ) 
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It follows that 



(47r) D 

0 . 


K 2 {H 2 aa') D ~ 2 r 2 (§) 


8 (D — 3)(D + 1) 


D 




( 112 ) 


4 RENORMALIZATION 


4.1 Counterterms 


To figure out what that appropriate BPHZ counterterms are, we are guided by two 
facts, the first of which is that we are calculating a 1PI graviton two-point function. 
This means our counterterms should contain two graviton fields. The second fact 
is the superficial degree of divergence of 4 in our primitive diagrams at one-loop 
order, which means that the counterterms should contain up to four derivatives. For 
quantum gravity at one-loop order, the necessary counterterms can be taken to be 
the Ricci curvature square R 2 and Weyl tensor square [33]. Note that 

we have obtained a perfectly de Sitter invariant representation of our graviton self¬ 
energy (152]) . thanks to the one-loop photon contribution being de Sitter invariant. As 
a result, there would be no additional counterterms. 

Moreover, it has been shown in section (13.2113.3ft that the nontrivial contributions 
of photons to the graviton self-energy come from the spin two part. Thus, the only 
counterterm we really need is the Weyl tensor square counterterm which is defined 
by the following Lagrangian: 


A£ 2 = ■ 


(113) 


Using the fact that the Weyl tensor with one upper index is conformally invariant, 



(114) 


We can rewrite the counterterm Lagrangian (11131) as 


A£ 2 = c 2 a D ~ 4 C a ^ 5 C a 01 s^ 

= c 2 a fl -V K /\^^(g/^d F )(Vd V ) + 0(A 3 ). (115) 


Note that we have used the linearized Weyl tensor operator C a ^ 1 s lll/ defined in section 
12.31 to expand the counterterm Lagrangian up to quadratic order in the graviton 
held. Then, one obtains the counterterm vertex by functionally differentiating i 
times counterterm action twice, and then setting the graviton held to zero, 


i5 2 AS 2 


2 c 2 ^rj aK r j ^ e rj^C a ^C' KXe / a x a D ~ 4 i 6 D {x - x') . (116) 


D-4-xD 


5h flu (x)5h pa (x') 


h=0 
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By recalling the definition (137)1 of J r ^ l/ P a one can see that the contribution of the Weyl 
tensor square counterterm to our graviton self-energy is 


i [^ASf ] (x; x') = 2c 2 K 2 T pvpa x 


a D A iS D (x 


x 


(117) 


4.2 Renormalizing the nonintegrability 

In the previous section, we derived our graviton self-energy in arbitrary D dimensions, 


- i [^E pa ] (x; x') = T pvpa F 2 (x; x') 


(118) 


where F 2 (x',x') is given by (lllip . One sees that our primitive result (I118p contains 
the factor of y 2 ~ D which diverges at coincidence (x M = x ,M ) with nonzero coefficients 
in D — 4 dimensions. Recall that to quantum correct the linearized Einstein field 
equation (JT)) we will integrate —i [ / "'E pCT ] over the four-dimensional measure; hence, 
(I118p is not integrable in those dimensions. However, only the terms that are at least 
as singular as y~ 2 need to be renormalized. Less singular terms are actually integrable 
and therefore can be evaluated in D = 4. To renormalize those nonintegrable terms, 
the general procedure is 

• First, convert them into functions of Ax 2 using the definition y = aa'H 2 Ax 2 . 

• Second, extract derivatives with respect to the unprimed coordinate x p , for 
example, 


1 d 2 1 

Ax 2 - 0-4 ~ 2 (D — 3)(.D — 4) Ax 2D ~ e ’ 

1 8 4 1 

Ax 2D - 2 ~ 4(D - 2) 2 (D - 3 )(D - 4) Ax 2D ” 6 ’ 


(119) 

( 120 ) 


in which the derivatives can be simply pulled out of the integration over primed 
coordinate x /p in the quantum corrected effective field equation (JT]). We keep 
doing this so as to get a less and less singular integrand until it becomes inte¬ 
grable in D = 4, but the price we pay for it is to get a manifest divergence, the 
factor of (D — 4) -1 . 

• Third, localize the explicit divergence of (D — 4)~ x in the previous step by 
adding zero using the identity 


d 2 


1 

Ax D ~ 2 


i4TT D / 2 

r (f -1) 




( 121 ) 
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This means that we can write 


d 2 


i47T D / 2 /i D 4 5 D (x — x') 


(D- A)Ax 2D ~ e T(f-l) (D — 4) 

d 2 


/i 


D—A 


(D — 4) Ax 2D - 6 Ax D - 2 


(122) 


where the factor of p with inverse length dimension is added for dimensional 
consistency. Then we expand the terms in the curly brackets away from D = 4, 




D—A 


(D — 4) in (/i 2 Air 2 ) 


Ax 2 


+ 0 


Ax 2D ~ e Ax D ~ 2 2 

So, we can write Eq. (11221) as 

<9 2 1 i47r°/ 2 jj, D ~ 4 8 D {x — x') 

(D - 4) Ax 2D ~ 6 ~ T (f - 1) (£) -4) 


(D- 4) 5 


(123) 


a 2 

' 2 


in (n 2 Ax 2 ) 
Ax 2 


+ 0(D — 4) . (124) 


Note that in Eq. (I124j) the logarithm term is integrable and the delta function 
piece is local though singular in D = 4. This is what we mean by localizing the 
divergence. 

• Finally, use the technique introduced above to break our result of graviton 
self-energy into finite and divergent parts. Then we shall set appropriate coun¬ 
terterm coefficients so as to cancel all the divergent parts, and obtain the fully 
renormalized result. 


Now, we are all set to renormalize the primitive result of — (11181) actually, 

the primitive result of F 2 (x]x') (1 111 1) . to be exact. Applying the identities (1119|) and 
(11241) in dm)) gives 


F 2 (x-,x') = 


hi 


6407T 4 


d 2 


ln(/i 2 Ax 2 ) 

Ax 2 


K 2 fi D ~ 4 T{ f) D(D - 2)i5 D (x - x') 


1287r? 


{D-A)(D-3) 2 (D + 1) 


(125) 


Substituting (11251) into (11181) and comparing it with (II171) . one sees that the local 
divergent part will be completely subtracted off by the counterterm contribution if 
we choose the arbitrary counterterm coefficient c 2 to be 


p D " 4 T(§) D(D- 2) 

2567rf (D ~ 4)(T> — 3) 2 (D + 1) 


(126) 
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Now, using this counterterm coefficient C 2 , we finally obtain the renormalized form 
of our graviton self-energy, 


i (x; x') = -i [ pu E p<T ] (x- x') - i ['“'ASf] (x; x') 


= T pvp<J hi (a)t5 4 (x - x') + F 2 , re (x; x') 
And the nonlocal part of our renormalized spin two structure function is 


F 2 tTe (x;x') = 


hi 


6407T 4 


d 2 


ln(/i 2 Aa; 2 ) 

Ax 2 


(127) 


(128) 


It is worth mentioning at this point that the spacetime dependence of equation (II 28f) 
allows us to easily reflect any derivatives on it, d' —> —d. This fact facilitates a dra¬ 
matic simplification of the tensor projector J^P u P a when acting on F 2 )Ve (Ax 2 ) in (1 127 [1 . 
which will be be useful in section mi To be specific, given the explicit expression for 
the projector F pupa in the Appendix of [32], one can see that in D = 4 dimensions it 
is reduced to the nice flat space form 

T pvpry ->■ i ^iP^lU 1 ' - ^W v W a ^j , IP 1 ' = d p d v - 7 f u d 2 . (129) 


Before closing this section, please note another fact, that with the following two 
identities, 


1 

Ax 2 
In (/i 2 Ax 2 ) 
Ax 2 


4 In (/.rAx 2 ) , 

(t 111 2 (/f 2 Ax 2 ) — 2 111 (/i 2 Ax 2 ) 


(130) 

(131) 


we can achieve a significant simplification by extracting another d’Alcmbertian from 

(H2HD 


F 2 , Ie (x-,x') = 


hi 


51207T 4 


d 4 


In 2 (/ 1 2 Ax 2 ) — 2 In (/i 2 Ax 2 ) 


(132) 


This form of the structure function (I132|) will come in handy when used in the 
Schwinger-Keldysh formalism. 


5 SOLVING THE EFFECTIVE FIELD EQUA¬ 
TION 

5.1 Schwinger-Keldysh formalism 

The in-out result of the graviton self-energy (I127|) seems to be ready for use in the 
linearized effective field equation (jT]) , but that would lead to two problems: 
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• Causality - The in-out graviton self-energy (j!27[) is nonzero for points x' 11 which 
lie in the future of x^, or at spacelike separation from it; and 


• Reality - The in-out graviton self-energy (11271) is not real. 


Note that there is nothing wrong with correcting the graviton propagator with the 
in-out graviton self-energy for asymptotic scattering amplitudes in flat space. But the 
causality and reality require that in the effective field equations, the appropriate 1PI 
two-point function should be that of the Schwinger-Keldysh (S-K) formalism. This 
technique provides a way of computing causal and real vacuum expectation values 
that is almost as simple as the Feynman diagrams which produce in-out matrix 
elements [31, S3136) 133 S3139, 40. @TJ . In our particular case, the photon propagator 
depends upon the y(x]x') which is a function of the Lorentz interval Ai 2 (x;x'). 
The four Schwinger-Keldysh propagators can be obtained by making the following 
replacements for the Lorentz interval [42]: 


Ax ++ — 

—* —*f 

x — x 

2 / \ 2 
— Urf-rfl-iS) , Ax 2 + _ = 

x—x 

2 ( \ 2 
- (r]-r]+i6J , 

III 

i 

CN 1 

< 

x—x' 

~ (\v~v'\+ity , Ax^ + = 

x—x' 

— (rj—rj'—iS'j . (133) 


To resolve the problem of reality and causality of the in-out formalism while used 
in an effective field equation, what we should really be using in (Q]) is the retarded 
graviton self-energy [^E^ 7 ] of the Schwinger-Keldysh formalism, 

[^E pa ] (x; x') [^Eg 7 ] (x; x') = ['"'E] (x; x') + [^Ef_] (x; x') . (134) 


It’s actually very simple to convert our in-out result (I127P into a S-K retarded one. 
We will get [ tlu E p ° + ] and ['“T)!] from (I127p by replacing the Lorentz interval Ax 2 
with Ax+ + and Ax+_. We will also drop the delta function terms in the -|— case 
and introduce an overall minus sign. Using (I132p we find that, 


r = k 2 ['“'E£’] ( 1 ) (z;x') 




ln(a) 

’80vr 2 


5 4 (x] x') + iF^^x; x ') 


with the S-K retarded structure function for the nonlocal part, 


iF 2 %(x] x') = 


id A 


5120vr 4 


In 2 (/i 2 Aa;+ + ) — In 2 (]i 2 Ax 2 + _) 


In ( fi 2 Ax 2 ++ ) — In [fi 2 Ax^_) 


(135) 


(136) 


Now, using the identities, 

In (/i 2 Ax^ + ) — In (/i 2 Ax 2 + _) 
In 2 {y 2 Ax 2 ++ ) - In 2 (/i 2 A4_) 


2 in6 (At] — Ax) , 

4 in6 (Ay — Ax) In y 2 (A if 



(137) 

(138) 
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we can achieve a manifestly real and causal form for iF^, 




-d A 

1280tt 3 


9 (A rj — Ax) In [/i 2 (A/? 2 — Ax 2 )] — 1 


(139) 


Here and henceforth we define conformal time intervals A r) = rj — rj', spatial intervals 
Ax = x — x', and its Euclidean norm Ax = ||Ax||. And 6 (...) denotes the Heaviside 
step function. 


5.2 Perturbative formulation 


To get the correction to force law of gravity, we need to solve the linearized effective 
held equation (ITJ) . However, we have only solved the lowest loop result for the graviton 
self-energy (I135j) . 

[^](x;x') = 0 + k 2 [^] {1 \x-x') + 0{k a ) . (140) 

Because the quantum corrections are only known to a finite order in k 2 , there is no 
alternative to making a similar expansion for the gravitational potentials, 

V(x) = h$ (x) + k 2 /iW(x) + 0(k a ) . (141) 

Substituting (I140j) and (I14ip into ([I]) leads to the zeroth and first order equations, 

= XM. (1«) 

V^Kh^ix) = S^(x) . (143) 

where we have the stress tensor density for the quantum corrections as 

S pv (x) = J d A x'[^Y, p £] { 1 \x-x')KhW(x') . (144) 


Note that we have regarded the matter source as zeroth order, unless its stress tensor 
includes loop corrections from the 1P1 one-point function. And the Lichnerowicz 
operator in de Sitter space is given by 


V' wpcT =-a 2 
2 


+Ha d 


— r) p,u i ) pa ) d 2 + rj pu d p d a + rf a d p d v — 2 

(rTrT - rf^Tf )v ) do - 2 ?f u 6 { 0 p d a) + 2S ( 0 p 7f Kp d u) ] + 3i/ 2 aV"W • (145) 


The force law of gravity is determined from the response of a point mass. So, consider 
a stationary point particle with mass M on a de Sitter background. We can write 
the matter action as 

^matter = ~M J dr - g ^(q{r))q^ {r)q u (r) 

= -M J dr|a(g(r))y/l - Kh 00 (g(r))| . (146) 
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To get from the first line to the second line in (11461) we have used the full metric (J5]) and 
have assumed the point particle has the world line q ^ = t5q. So, the corresponding 
linearized stress tensor density takes the form 


T^(x) = -k 


^matter [h\ 


5h flu (x) 


h =0 


— 87tGMSq()q x<u5 3 (x) . 


(147) 


With this matter source, the solution to the zeroth order equation 111 Id is given by 

1251 

m + lev) 2 4¥r ■ (148) 

a i \ x 11 

Then, integrating 1 up against this tree order solution hp,? (j!48j) simply gives 

the one-loop source to the first order equation (1143|) which is to be solved. 


5.3 The one-loop source term 

The quantum source il I 1 Hi consists of both local and nonlocal parts, 

^(x)=5C: al (,) + 5l local (x). (149) 


The local part can be very simply expressed in terms of the linearized Weyl tensor 
(juavp^x) formed from nhfj [32] . 


= d a dp 

The nonlocal part of the source, 


ln(a) 

407r 2 


(JIMXV 0 | 


XI 


dL„c.,M = I dlx ' 'dXfo 1 ') Kft ISM, 


(150) 


(151) 


requires two steps: performing the integral and acting the derivatives. 


• Fundamental integral 


Substituting fll39[) and (I148jl into (11511) implies that there is a fundamental integral 
to do, 


F(V, r) 


= / dV< 


6 (Ar] — Ax) 


In 


/r 2 (Arf — Ax 2 ) 


- 1 


(152) 


After changing the spatial integration variable, x! —> f = Ax, letting x = ||x||, and 
performing the angular integration under spherical polar coordinates, one sees that 
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the integral (1152)1 becomes 


2vr 

x 



(x + r 


x — r |) rdr 


x f 6(Arj — r) 

•J Vi 


In 


A 


(Ar/ 2 




(-ddrj')drj' , 


(153) 


where rji = —H^ 1 denotes the initial time. Now, let us first tackle the temporal 
integral. After using the Heaviside step function to fix the temporal integration 
limits, we find 



- r) 

In 

n~ r 

In 

'Vi 



ji 2 (A rj 2 — r 2 ) 
ji 2 (A rj 2 — r 2 ) 


- 1 


(—ddrj')drj' 

'—ddrj')drj' 


= 2dd (Arji + r}i)r ln(2 fir) + dd(2Arji + 3 rji - r)(Arji 
~H(Ar)i + 2r]i + r)(Arji + r) ln[/_t(A^ + r)] 

~^H(Ar)i + 2rji - r)(Arji - r) In [//(A rji - r)] 


r) 


(154) 


where we have defined A rji = rj — rji. Note that the result (1154ft will force the 
integration limit of the radial integral in (1153ft to be 0 < r < A rji. Next step is to 
substitute (1154)1 into (11531) . and we are left with a purely radial integration which 
is tedious but straightforward to do. Here, we simply give the full result of the 
fundamental integral (11521) . 


/ \ 2tt dd I / a \ 4 / 

F(rj,x) = x < -( Arji + rji)x ln(2 jix) 

5(119A77j + AAQiji)Aijfx + 1250(A?7j + rji)x A — 10(131Ar/j + 220rji)ArjiX 3 — 261x 5 


600 


20 L 


(A rji + 10 rji)x — 2(A rji + 5rji)Arji + 3x 2 (Arji + x) 3 ln[/i(A?/j + x)] 


20 


(Arji + 10rji)x + 2(Arj t + d>iji)Arji - 3x 2 (Arji ~ xf In [ji(Arji - x)] \ . (155) 


• Acting derivatives 


At this stage, we can rewrite the nontrivial part of our nonlocal source term (I15ip as 


^nonlocal (V, *) 


GM 

6407T 3 


(2 5 0 p 5 0 a + rj p(T )T^d 4 F( V ,x) . 


(156) 
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Next is to act two d’Alcmbertians on F(r),x ) using the following differentiation for¬ 
mula: 

<9 4 F(?y, x) = F (4 ’ 0) - 2F (2 ’ 2) + F (0 ’ 4) - -F (2 ’ 4) + -F (0 ’ 3) , (157) 

xx 

where we have adopted a shorthand notation of mixed partial differentiations, 


ain+n 

Fb "- n) ^W^ F M- 

Then, it follows that two d’Alembertians on (j 15 5ft gives 

GOq, x) = d 4 F(‘r], x) = IQixFl 1 —-In (2 nx) 

L x 


(158) 


(159) 


Now the remaining work is acting the operator (2h p <5° + r/ p(J ) on G(rj, x). Note 

that the tensor projector F 7 ^' 7 is traceless by design; therefore, what we really care 
about is the result of acting (2<5 p <5°) F p,jp(J , 

28 ( XF llupa = d p d u d 2 0 - -d' J d v d 2 -28^d u) d 2 d 0 +-rj" v d 2 d 2 + (s p S" + -r/'A 8 4 . (160) 


Similarly, we can also derive a formula for this differential operation, 

2S°XF pupa G{r), x) = d“d v (g {2 ’ 0) - iG (0 ’ 2) - —G (0,1 A 
' \ 3 3 x J 

+5ffi f--G (1 ’ 2) - — G^A + 8 P 5 U 0 ( -G (0 ’ 4) + — G (0 ’ 3) 

V 3 3 oc J V 3 3 oc 

+Tj lw f--G (2 ’ 2) + -G (0 ’ 4) - — G (2il) + — G (0 ’ 3) ^ . 

' V 3 3 3 x 3 x J 


(161) 


Substituting the result of G(rj,x) (11591) into (11611) allows us to finally obtain the 
nonlocal part of our one-loop source term, 


5 nonlocal(^) = 


GM f— n— v (lh] 


407T 2 | 


8T8 


3x 3 ' Mo 


4Hr] 


x u 




(162) 

Substituting (11501) and (11621) into (j!43j) . and ignoring the delta function contributions, 
gives the following nontrivial components of our first order equation: 


GM 


"-M) = 1} . 

= a ‘ 9 ’{-T^ + l^} + sv {^}- 


(163) 

(164) 

(165) 
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5.4 Quantum corrected potentials 

We want to express the one-loop Newtonian potentials in Schwarzschild coordinates 
so that with a convenient gauge choice their nonzero components take the form 

«/*oo (v, x) = h (rj, x) , rJi^ ( 77 , x) = 0 , Kh\f(r},x) = f 2 (r],x)5 ij . (166) 

Substituting our ansatz (11661) into (1163!) - (11651) gives rise to the following equations 


GM GMH 2 

607r 2 a 3 a: 3 207r 2 aa; 

GM 


607 r 2 a 3 x 3 
GM 

107r 2 a 3 x 5 

GM 

20ir 2 a 3 x 5 


fi — f2 , 

r,if 0) ~ h , 

/f 2) + -M + -fi m 

x r] 



-Ji 


( 0 , 2 ) 


+ -/i (0,1) - -/r 

z X 1) 

_ 1 ,( 0 , 2 ) _ 1 ,( 0 , 1 ) _ 2 ( 1 , 0 ) ( 2 , 0 ) 

2 J 2 J 2 J 2 t/ 2 

X Tj 


- (i ' o) +4/. 

T] z 


(167) 

(168) 
(169) 


(170) 


We should remind readers of the notation for mixed partial differentiation (I158j) . Now 
it seems that we have four equations (I167|) - (1170|) but only two field variables, f\ and 
/ 2 - Actually, these equations are not completely independent. We can pick the easier 
ones to solve; then, the remaining equations shall be automatically satisfied by the 
virtue of stress energy tensor conservation and the Bianchi identity. 

Equations (1 1671) and (I168j) are much simpler to start with, so we first eliminate 
fi using these two and obtain a first order partial differential equation for / 2 , 


^ 0/2 + Haf 2 


GMH GMH 3 

307r 2 a 2 ic 3 207r 2 a 


whose solution is 

ft A- GM GMH 2 ln(a) C{x) 
' 2 \ r h x ) 607r 2 a 3 ic 3 207r 2 aa: ^ a 


(171) 


(172) 


Note that C(x) is a integration “constant,” but it may still depend on the Euclidean 
norm x. To fix this constant, we can substitute (I167|) and (II 72[) into (I169j) and obtain 
a second order ordinary differential equation for C{x), 


C"(x) + ~C\x) 
x 


GMH 2 

207r 2 a: 3 


(173) 


which can be regarded as first order in C{x) and hence easily solved. Thus, we find 


C{x) 


GMH 2 
207 t' 2 x 


1 + In ( Hx ) 


+ — + B, 
x 


(174) 
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where A and B are really constants. To fix them we shall, in principle, make use of 
(11701) . Well, as we have mentioned previously, whichever equation is left should be 
satisfied automatically. And indeed, with the solutions given by (I167f) and (11721) and 
(11741) . one can show that (11701) is solved no matter what A and B are. Therefore, we 
can simply set A — B — 0, and we obtain the formal solutions for f\ and / 2 , 


h{rhx) 


GM GMH 2 In (aHx) 

307 f 2 a 3 x 3 207r 2 aa; 

GM GMH 2 [l+\n(aHx)) 

607r 2 a 3 x 3 207r 2 ax 


(175) 

(176) 


Combining the tree order solutions (I148p gives the final results of our one-loop New¬ 
tonian potentials, 
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2 



K, 


K 


^ h 0i (v, x) 
( T ], x) 


K 


60n 2 (ax) 2 

0(k 4 ) , 


+dAhM£) +0(K -‘ ) U kC’ 6.2) 


J _ a _ : 

|^1207r 2 (aa:)^ 


407r 2 


+ dA!i±!A2Ci+o (K ‘ ) U 


407T 2 


(177) 

(178) 

(179) 


6 DISCUSSION 


We have computed the one photon loop contribution to the graviton self-energy on de 
Sitter space, and we have used it to derive quantum corrections to the gravitational 
potentials of a stationary point mass. Our results (II77[) ~ (II79|) feature two sorts of 
terms which can best be understood by expressing the full Newtonian potential as 
its classical value times a series of corrections, with the factors of h and c restored, 


<h(?7, x ) 


2 GM 
ax 


1 + 


4 hG 

157rc 3 (ax) 2 


+ 


2 hGH 2 
57rc 5 



+ 1 


+ 0(G 2 


(180) 


The fractional correction 4frG/[157rc 3 (a:c) 2 ] is just the de Sitter version of the flat 
space background correction first found in 1970 by Radkowski [33] and later con¬ 
firmed by Capper, Duff and Halpern [ 33 ] • Getting this term right checks an impor¬ 
tant correspondence limit; however, it can only be significant at physical distances of 
ax ~ 10 35 m. This is a small distance even with respect to the scales 

prevailing during primordial inflation. The more interesting part of (I180|) is the frac¬ 
tional correction 2hGH 2 /[5irc 5 ] x In (aHx/c), which vanishes in the flat space limit. 
This represents a small logarithmic running of the Newton constant. One interest¬ 
ing consequence is that an observer carried away by the Hubble flow at constant x 
sees the one loop correction to the gravitational potential grow with respect to its 
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classical value, although classical and quantum potentials both fall off exponentially 
in co-moving time. Hence a very long phase of inflation leads to a breakdown of 
perturbation theory. 

The logarithmic running correction in (11801) might seem surprising because it 
is sourced by the gravitational vacuum polarization induced by photons, which is 
simply related to its value on flat space by the conformal invariance of classical 
electromagnetism. Indeed, the source terms we found in Eqs (I163l) - (ll65j) are almost 
identical to those of flat space. The sources on the right-hand side of (I163|) and 
(11651) are just conformal rescalings of the flat space result, and the new, “de Sitter” 
source on the right-hand side of (11641) is the bare minimum required by conservation in 
view of the time dependent conformal rescaling. The reason this conformally rescaled 
source can do more than “de Sitterize” the response of the flat background is that the 
left-hand side of the equation is different because gravity is not conformally invariant. 
Gravity in a de Sitter background responds far differently to sources than gravity does 
in a flat background. 

It is worthwhile working out the response seen by a static observer whose back¬ 
ground geometry is 


ds 2 = - (l - H 2 r 2 ) dt 2 + 


dr 2 


+ r 2 dil 2 


(1 - H 2 r 2 ) 


The transformation to these coordinates is 


-1 „-Ht 


-H~ 1 e 


V = 


x = 


y/1 -H 2 r 2 


/y'i' g Ht 


\J\-H 2 r 2 


dri = — 
a 


dx l = - 
a 


dt 


Hrdr 

1—H 2 r 2 


-Hr l dt + dr 1 + 


H 2 r l rdr 
1 -H 2 r 2 


(181) 


(182) 

(183) 


The potentials come from using transformations (I182j) - (11831) on the full metric g^ u = 
a 2 [r]^ + Khpvl, 


7 2 _ 2 

(ZfSfuji — CL 


-dr] 2 + dx-dx + Khooa 2 di] 2 + 2nhoia 2 dr]dx l + nh^a 2 dx l dx 2 , (184) 

dr 2 


1 -H z r 


2„2 


dt + 


1 — H 2 r 2 


+ r 2 dfl 2 -\rXttdt- J + 2xoidtdr l + Xijdr l dr J . (185) 


The final results are complicated by the fact that achieving the de Sitter-Schwarzschild 
form would require order GM modifications of the transformations (I182|) and (I183p . 
However, one important feature of the result is that all of the factors of ax become the 
time independent static coordinate radius r. Hence an observer at constant physical 
distance from the source detects only a temporally constant, logarithmic enhancement 
of the potential. 
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